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摘要 : 大气 中 水 蒸气 会 影响 气 溶胶 颗粒 的 生长 演化 动力 学 特性 , 从 而 改变 其 环境 效应 。 当 前 ， 对 包含 有 机 组 分 气 溶胶 颗粒 与 

分 子 动力 学 模型 ， 分 

重点 分 析 了 各 影响 因素 

作用 下 团 徐 结构 和 气 粒 界面 特性 。 结 果 表 明 ， 温 度 升 高 ，CaHi01-H20 团 艇 相继 出 现 分 层 结 构 和 混合 结构 。 加 入 40 个 卤化 钠 分 

子 后 ， 水 分 子 对 初始 团 簇 溶解 程度 的 排序 为 40 NaF <40 NaCl < 0 Nax <40 NaI， 且 差异 随 着 温度 的 升 高 而 变 大 。 高 温 ( 企 300 
K) 时， 各 种 团 簇 的 表面 均 包 含 CH0, 分 子 和 HO 分 子 ， 且 CH0, 分 子 的 疏水 基 朝 向 气相 一 侧 。 


水 相互 作用 机 理 的 认识 还 很 局 限 。 本 文 建立 了 纯 有 机 气 溶胶 和 有 机 -无 机 混合 气 溶胶 颗粒 水 溶 特性 
别 模拟 了 不 同 温度 条 件 下 不 同 讽 化 钠 盐 (NaI/NaCl/NaF) 对 丙 二 酸 (CH0) 团 复 水 溶 特性 的 影响 情况 。 
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The Impact Mechanism of Sodium Halide on the Water Soluble Characteristics 


of Malonic Acid Aerosol Particles 


ZHANG Chao, WANG Yue-She, LI Chen-Pei, WANG Hai-Jun 
(State Key Laboratory of Multiphase Flow in Power Engineering, Xi'an Jiaotong University, Xi’an 710049, China) 


Abstract: Water vapor in the atmosphere can affect the growth dynamics of aerosol particles, thus changing their environmental 
effects. However, the understanding of the mechanism of the interaction between water and aerosol particles containing organic 
components is still limited. In this study, molecular dynamics model is established to study the water soluble characteristics of pure 
malonic acid (C3H4O,) aerosol and malonic acid-sodium halide mixed aerosol particles at the condition of different temperatures and 
different sodium halide types (NaI/NaCl/NaF). The key point is the effect of different factors on the structural and the interfacial 
properties of particle. The results show that the demixed structure and the mixed structure can be observed for malonic acid-water 
aerosol with the increase of temperatue. After adding 40 sodium halide, the order of the dissolution degree of the initial particle is 40 
NaF < 40 NaCl < 0 NaX < 40 Nal, and the difference increases with the increase of temperature. At higher temperatue (T = 300 K), 
the interface of various aggregates is conposed of C3H4O4 and H2O molecules, and the hydrophobic group of C3HsO4 molecules faces 


the gas phase. 
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0 引言 


大 气 溶胶 颗粒 物 会 影响 全 球 气候 状况 、 大 气 辐 对 气 溶胶 颗 粒 的 光 散 射 特性 ， 云 寿命 和 颗粒 沉降 率 
射 强 迫 和 区 域 性 空气 质量 。 环 境 中 水 蒸气 会 对 气 溶 的 影响 ““。 此 外 ， 有 机 物 在 气 溶胶 颗粒 的 组 成 成 


胶 颗 粒 物 生长 演化 过 程 产生 重要 作用 ， 主 要 表现 为 分 中 占有 很 大 比例 ”“， 其 与 水 燕 气相 互 作用 时 的 动 
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力学 演化 过 程 一 直 是 气 溶 胶 研 究 领 域 的 热点 和 难 
点 。 目 前 已 有 的 研究 主要 采用 实验 方法 和 基于 热力 
学 经 典 理论 建 模 的 方法 。 这 两 种 方法 在 研究 气 溶胶 
和 水 相互 作用 过 程 中 分 子 层面 机 理 问 题 上 有 很 大 局 
限 ， 因 此 需要 借助 分 子 模拟 手段 。 

当前 ， 基 于 经 典 力学 理论 的 分 子 动力 学 方法 已 
经 逐渐 被 应 用 到 简单 成 分 气 溶胶 颗粒 吸湿 特性 动力 
学 过 程 的 研究 中 。 如 ，Zhang 等 人 “采用 MD 方法 
研究 了 均 相 成 核 和 离子 诱导 水 蒸气 成 核 过 程 中 新 生 
成 团 复 演化 过 程 的 动力 学 特性 。Bahadur 和 
Russell 采用 MD 方法 研究 了 NaCl 颗粒 的 潮解 特 
性 并 分 析 了 吸湿 过 程 中 水 分 子 吸附 系数 的 变化 过 
程 。 Akbarzadeh 等 人 "采用 MD 方法 研究 了 NHC1 
和 NENO: 颗粒 潮解 特性 的 差异 。Picaud 课题 组 ”“ 
基于 MD 方法 分 别 研究 了 丙 二 酸 、 甲 酸 、 乙 酸 和 丙 酸 
气 溶胶 颗粒 与 水 分 子 相互 作用 下 二 元 气 溶 胶 颗 粒 结 
构 的 演化 过 程 ， 据 此 分 析 了 不 同 种 类 二 元 颗粒 的 吸 
湿 特 性 。 在 将 模拟 结果 与 Schill 的 实验 工作 "“ 相 
对 比 后 得 出 结论 : 有 机 气 溶胶 颗粒 的 吸湿 特性 不 仅 
取决 于 有 机 物 中 的 氧 碳 比 ， 而 且 还 与 有 机 物 中 羧基 
数量 和 有 机 分 子 的 结构 相关 。Chakraborty 课题 组 
一 呈 先后 研究 了 一 元 羧 酸 和 二 元 羧 酸 包 庄 水 滴 气 溶 
胶 颗 粒 的 结构 和 界面 特性 。 然 而 ， 当 前 对 有 机 物 和 
无 机 盐 混合 气 溶胶 颗粒 物 和 水 蒸气 相互 作用 过 程 的 
研究 主要 采用 实验 方法 ， 对 其 分 子 层面 动力 学 机 理 
的 研究 还 很 缺乏 ” 。 

本 文 建立 了 丙 二 酸 (CHO) -水 (H0) 二 元 体系 
All A) = (C3H.0.) - Be 14,44] (NaI/NaCl/NaF)- 水 (上 0) 
三 元 体系 的 分 子 动力 学 模型 ， 模 拟 了 不 同 温度 条 件 
下 ， 不 同 种 类 亢 化 钠 盐 的 存在 对 有 机 气 溶胶 颗粒 和 
水 分 子 相互 作用 过 程 的 影响 ， 对 比分 析 了 颗粒 物 的 
结构 和 界面 演化 特性 。 


1， 分 于 动力 学 模拟 
本 文 模拟 体系 中 包含 三 种 分 子 ， 其 中 水 分 子 采 
用 SPC/E 模型 , 丙 二 酸 分 子 采用 OPLSAA 力 场 中 的 势 


能 参数 ”“， 卤 化 钠 中 的 离子 采用 文献 ”中 势能 参 
数 。 模 拟 体系 中 势 函 数 表达 式 为 : 
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其 中 ， 公 式 右边 第 一 项 为 非 键 接 势能 项 ， 包 括 
Lennard-Jones 12-6 势能 和 静电 势能 ,对 于 同一 个 
分 子 中 1-4 原子 相互 作用 而 言 , A0. 5, 否则 , EFL. 
后 三 项 为 分 子 内 相互 作用 势能 项 , 分 别 为 键 伸缩 项 ， 
键 角 弯 曲 项 和 双 面 角 扭曲 项 。 范 德 华 作用 力 和 库伦 
作用 力 截断 半径 为 1. 4 nm， 长 程 库仑 力 采 用 PPPM 
算法 。 模拟 采用 周期 性 边界 条 件 , 时 间 步 长 为 1 fs. 
本 文中 的 模拟 过 程 均 采 用 NVT 系 综 ， 温 度 控制 采用 
Nose-Hoover 热 浴 。 
首先 ， 建 立 三 个 分 别 包含 120 个 丙 二 酸 分 子 
(CsH:04) 和 40 个 不 同 卤化 钠 分 子 (Nal/NaCl/NaP) 的 
团 簇 来 分 别 不 同 有 机 -无 机 混合 气 溶胶 颗粒 ,将 其 放 
置 于 边 长 为 8 nm 的 立方 体 盒子 内 。 经 过 2 ns 的 平 
衡 后 将 团 簇 转 移 到 一 个 新 立方 体 盒子 ( 边 长 8 nm) 的 
心 ， 同 时 加 入 360 个 均匀 分 布 的 水 分 子 来 模拟 一 
定 湿度 条 件 。 体 系 在 100 温度 条 件 下 持续 5 ns. 
在 此 过 程 中 发 现 所 有 模拟 体系 中 的 水 分 子 都 被 吸附 
到 初始 团 秘 上。 此 时 ， 将 团 徐 再 次 放 入 一 个 新 的 立 
方 体 盒子 中 ( 边 长 22 nm), 并 在 不 同 温度 条 件 〈100， 
140, 180, 220, 260 和 300 K) 下 持续 1 ns。 本 文 
所 采用 的 一 些 温度 值 低 于 大 气 对 流 层 的 温度 ， 这 是 
因为 在 分 子 动力 学 模拟 中 ， 团 簇 结构 和 界面 的 变化 
主要 依赖 于 所 采用 的 分 子 势能 参数 ， 由 于 分 子 势能 
参数 的 精确 性 不 能 保证 ， 因 此 在 满足 团 簇 分 子 不 燕 
发 的 前 提 下 将 模拟 体系 控制 在 不 同 的 温度 ， 以 获得 
团 簇 可 能 出 现 的 所 有 结构 和 界面 特性 。 最 后 ， 对 各 
模拟 体系 进行 1 ns 的 数据 采集 。 

作为 对 比 ， 本 文 进行 了 纯 丙 二 酸 团 徐 (120 
CH0) 的 吸湿 性 模拟 ， 模 拟 过 程 同 上 。 为 了 分 析 不 
同 团 簇 的 结构 和 界面 特性 ， 在 分 子 模拟 结果 的 基础 
上 ， 本 文 分 别 分 析 了 各 体系 中 团 秘 平衡 快照 、 团 簇 
尺寸 分 布 函 数 和 团 簇 径 向 分 布 函 数 。 分 析 过 程 中 ， 
采用 Stillinger “标准 判断 两 个 分 子 是 否 相 连 ,不 
同 种 类 分 子 间距 离 标准 分 别 为 : 两 个 两 二 酸 分 子 
(满足 其 中 任 一 条 件 即 为 相连 ) : (0 (C=0) -H (0-H) ) 
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< 0.225 nm, r(0(C=0)-H(C-H)) < 0.325 nm, 
r(0(C-0)-H(O-H)) 和 0.225 nm, r(0(C-0)-H(C-H)) 
<0. 325 nm。 丙 二 酸 分 子 和 离子 : r(Na’-0(C=0)) < 
0. 32 nm, r(Cl -H(O-H)) 和 0.28nm, r(F-H(0-H)) 
< 0.19 nm, r(T-H(@-H)) < 0.325 nm. BAI 
离子 : r(Na-Cl) < 0.33 nm, r(Na-F) < 0.28 
nm, r(Na-Il) < 0.41 nm， 这 些 临 界 值 为 相应 径 
向 分 布 函数 第 一 个 波 谷 对 应 的 距离 值 。 


2. 模拟 结果 与 分 析 


2.1 团 簇 结构 图 

不 同 模 拟 体 系 中 团 簇 的 结构 可 以 从 团 簇 模拟 快 
照 中 获取 定性 结果 。 本 文 分 别 追 踪 了 不 同 温度 条 件 
下 各 种 类 型 团 复 的 结构 变化 情况 。 图 1 中 显示 了 两 
种 温度 (140 K 和 300 K) 条 件 下 纯 CHO, AR 
/CsHi0.+Nal/NaC1/NaF 团 徐 和 水 相互 作用 后 的 结构 
快照 。 从 图 中 可 以 看 出 ， 对 于 二 元 CH40s+H20 AA 
而 言 ， 随 着 温度 的 升 高 ， 团 簇 会 相继 出 现 分 层 结构 
(水 分 子 以 团 复 形式 包 庄 在 CHO ARID) 和 混合 结 
构 (水 分 子 渗入 CH0, ANKE) o 。 此 二 元 C3H,0.+H,0 
团 徐 的 变化 情况 和 文献 ”中 模拟 结果 相同 ， 定 性 证 
明了 本 文采 用 的 模拟 方法 的 正确 性 。 

当 加 入 元 化 钠 后 ， 随 着 温度 升 高 ， 三 元 
CaHs0stNaI/NaCl/NaF+Hs0 B fe 2 HH BS = 7G 
CsHs04+Hz0 AeA AG. 其中, 加 入 Nal 和 NaCl 
后 团 复 会 出 现 分 层 结构 和 混合 结构 , 而 加 入 NaF 后 ， 
团 复 只 出 现 分 层 结构 。 


= 


3 0 Nax 40 Nal 3 
=| a k ii, 
Ž 40 NaCl 40 NaF 
z Re v ay 
140 K 300 K 
图 1. 不 同 温度 (140, 300 K) 条 件 下 CHO+Hz0 团 徐 和 
CsHiO;+Nal/NaC1/NaF+H.0 团 簇 的 模拟 快照 。 


Fig.l Snapshots of aggregate containing C3H4O4 and H,O and 
aggregates containing C3H4O., NaI/NaCl/NaF and HO under 
the condition of diffrerent temeperatures (140/300 K). 


2.2 ARR AHR 
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不 同 温度 条 件 下 二 元 CHO:+H0 HiME 
CsH:0.+Nal/NaC1/NaF+H.0 团 艇 结构 变化 情况 的 定量 
描述 可 以 通过 计算 模拟 体系 中 CHO A RE BK 
CsHi0.+Nal/NaC1/NaF 团 簇 的 尺寸 分 布 函数 获取 。 与 
图 1 相对 应 ， 图 2 列 出 了 两 种 温度 (140K 和 300 K) 
条 件 下 不 同类 型 团 簇 的 尺寸 分 布 函 数 (a-qd) 。 同 时 图 
2 也 给 出 了 温度 为 300 K 时 不 同 种 类 团 艇 尺寸 分 布 
函数 第 二 个 峰值 的 位 置 (e) 。 为 了 便于 比较 ， 在 图 
2(e) 中 将 二 元 体系 中 CHO 团 艇 的 所 有 尺寸 在 原来 
的 基础 上 分 别 乘 以 5/3。 


1.0 
0.8 
0.6 | 
S | 
ù 0.4 00 0 60 120 0 A 
0.2 
0.0 1 s a 
0 60 120 0 60 120 
n 
T=300 K 
T=140K T=300 K 


A 
009 0 60120180 0 60120180 4 
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(e) 
图 2 不 同 温度 (140 K, 300 K) Ae FRA BERRIA BR 


BL P(n), (a) CH0: ABE (b) CHOtNar AR, (c) 
CHO+NaCl HN, (d) CH:0O+NaF (图 中 插图 为 局 部 放大 图 ) , 
(e) 7300 K， 各 体系 中 团 复 尺寸 分 布 函数 的 第 二 个 峰值 
Fig.2 Cluster size distribution P(n) of various cluster, (a) 
C3H,O, cluster, (b) C;H4,O,+Nal cluster, (c) C3H,O4+NaCl 
cluster, (d) C3H,O,+NaF cluster (Insets of the figures show the 
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corresponding P(n) in different scale) , (e) at 7=300 K, the second 
peak of P(n) in each system. 


对 于 二 元 CaH404+Hs0 体系 中 的 CHO. ABET A. 
不 同 温度 条 件 (140 K 和 300 多 下 的 尺寸 分 布 函数 
都 存在 两 个 峰值 (n 守 1，120 处 ) ， 对 应 体系 中 主要 
出 现 的 两 种 团 簇 尺 寸 . 其 中 尺寸 近似 为 1 的 团 簇 对 
应 体系 中 略 超出 本 文 定义 距离 判断 标准 的 分 子 。 温 
度 升 高 ， 第 一 个 峰值 增 大 ， 第 二 个 峰值 减 小 ， 表 明 
水 分 子 渗入 初始 CHO 团 簇 ， 团 簇 被 部 分 溶解 。 这 
种 现象 的 发 生 是 由 于 分 子 热 运 动 的 加 剧 导 致 丙 二 酸 
团 簇 内 部 产生 了 大 量 空 穴 ， 致 使 水 分 子 渗入 其 中 。 
加 入 商 化 钠 盐 后 ， 两 种 温度 条 件 下 CHi04+NaX 团 簇 
尺寸 分 布 函数 的 波峰 数目 及 变化 趋势 与 二 元 体系 中 
CH0, 团 簇 基 本 相同 。 然 而 ， 由 图 2(e) 可 以 看 出 ， 
四 种 团 秘 分 布 函 数 第 二 个 峰值 所 对 应 的 尺寸 值 排 序 
为 40 Nal < 0 NaX < 40 NaCl < 40 NaF, 这 种 现象 
Be eA CAN ER A SRK TGR CH0, 团 簇 
EVA APE, m ALS EE AY OE ES La 
盐 的 种 类 有 关 。 此 结果 定量 验证 了 上 文中 团 簇 模拟 
快照 分 析 结 论 。 

为 了 更 加 直观 体现 不 同 因素 对 模拟 体系 中 初始 
CsHi01 BIFE CsH40s+Nal/NaC1l/NaF 团 簇 被 水 分 子 溶 
解 程度 的 影响 情况 ， 本 文 计算 了 不 同 模拟 体系 中 最 
KAERT M 相对 于 初始 尺寸 (WW) 的 变化 率 CAM) 
随 温 度 的 变化 情况 ， 结 果 如 图 3 所 示 。 


1.00 


—e— 0 NaX 
—@— 40 Nal | 
一 全 一 40 NaCl | 
—v¥— 40 NaF | 


0.85 


200 240 . 280 l 320 
T/K 
图 3 模拟 体系 中 最 大 CsHsO04 才 JEER CsHi0.tNaI /NaC1/NaF AiR 
的 变化 率 (WN) 随 温度 的 变化 情况 


Fig.3 Change rates (M/N) of the largest C3H4O4/ 


L 1 
120 160 


C3H4O4-NaLNaCLlINaF aggregates in different systems, as a 


function of the temperature 
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HE 3 可 以 看 出 ， 低 温 时 (7<180 K), MAT 
拟 体系 中 初始 团 簇 被 水 溶解 的 程度 近似 为 零 , 随 着 
温度 升 高 ， 所 有 模拟 体系 中 的 初始 团 簇 尺寸 减 小 ; 
相同 温度 条 件 下 团 簇 尺寸 的 变化 率 (AM) 排序 为 40 
NaF > 40 NaCl > 0 NaX > 40 Nal. XXi HH pakt 
盐 的 加 入 会 改变 初始 团 秘 的 稳定 性 ， 此 外 ， 枫 化 钠 
盐 对 初始 丙 二 酸 团 簇 稳定 性 影响 程度 的 差异 随 着 温 
度 的 升 高 而 增 大 。 


T 


2. 3 团 簇 径 向 分 布 函 数 

湿度 条 件 下 ， 团 复 中 不 同 组 分 的 分 布 情况 可 由 
相应 分 子 中 代表 性 原子 的 径 向 分 布 函数 获得 。 由 于 
组 分 径 向 分 布 函 数 的 计算 要 求 团 簇 为 近似 球形 ， 因 
此 本 文 只 计算 了 7180 K 和 300 K 条件 下 
CaHs0s+Hz0 fe A CsaHi0s+Nal/NaCl/NaF+Hz0 团 簇 中 
CHO 分子、E0 分 子 、Na 和 XX 的 径 向 分 布 函数 。 其 
中 ，0(C=0) 和 H(C-H) 代表 CHs04 分 子 ，0 (Hz0) 代 表 
Hz0 分 子 。 结 果 如 图 4 所 示 。 
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4 不 同 温度 (180 K, 300 K) 条 件 下 ， 二 元 CsH.0.-H.0 BGK 
a) 和 三 元 CsH:0.+Nal /NaC1/NaF-H.0 Hf (b, c~ d) 中 不 同 原 
子 的 径 向 分 布 函数 
Fig.4 Radial density profiles of different atoms in the binary 
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C3H4O4-H2O aggregate (a) and the ternary 
C;H40,+Nal/NaCl/NaF-H,0O aggregate (b/c/d), under the 
condition of different temperatures (180 K, 300 K) 


图 4 可 以 看 出 ， 对 于 所 有 模拟 体系 ， 商 化 钠 
盐 离 子 位 于 团 艇 内 部 。 较 低温 度 时 ( 产 180 kK) ， 水 分 
FAREYLE CH0 HEM CsHs0s+tNal/NaCl/NaF 
外 ， 随 着 温度 升 高 (六 300 K), KAT E ASEAN BH 
散 ， 但 相同 温度 条 件 下 ， 不 同 团 簇 中 水 分 子 扩散 程 
度 不 同 ， 可 以 明显 地 看 出 ， 包含 NaCl/NaF 的 丙 二 
酸 团 簇 中 水 分 子 的 扩散 程度 小 于 纯 丙 二 酸 团 徐 和 包 
S Nal 的 丙 二 酸 团 艇 。 且 在 较 高 温度 ( 堪 300 K) 时 ， 
对 于 二 元 CaHs0s+H20 团 徐 和 三 元 C:H0,+H:0+NaX 团 
fe “NF AAR FAK Fak, AAS 
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酸 分 子 中 的 疏水 基 (-CH) 朝向 气相 一 侧 ， 亲 水 基 
(—COOH) 朝向 颗粒 相 内 部 。 

气 粒 界面 不 同 组 分 的 面积 比 以 及 有 机 物 分 子 的 
形态 会 影响 气相 一 侧 水 花 汽 分 子 在 团 徐 表面 的 粘 附 
系数 "”"， 从 而 影响 该 团 簇 吸 湿 生 长 特性 。 通 过 本 文 
的 模拟 分 析 可 以 获得 纯 丙 二 酸 以 及 丙 二 酸 -卤化 钢 
混合 气 溶胶 颗粒 表面 不 同 成 分 相对 含量 的 变化 规律 
以 及 表面 有 机 物 组 分 的 分 子 形态 ， 这 将 为 后 期 气 游 
胶 颗 粒 表面 水 分 子 吸附 系数 的 定量 描述 黄 定 基础 。 


3， 结 论 


本 文采 用 分 子 动力 学 方法 模拟 了 纯 丙 二 酸 有 机 
气 溶胶 颗粒 和 丙 二 酸 -卤化 钠 混合 气 游 胶 颗 粒 的 水 
溶 特性 ， 通 过 模拟 快照 的 分 析 以 及 团 簇 尺 寸 分 布 、 
最 大 团 复 尺寸 变化 率 和 团 复 径 向 分 布 函数 的 计算 ， 
分 别 从 定性 和 定量 两 个 角度 研究 了 温度 和 讽 化 钠 盐 
种 类 (Nal/NaCl/NaF) 对 丙 二 酸 (C0i) 气 溶胶 颗粒 
水 溶 特性 的 影响 规律 和 微观 机 理 。 主 要 研究 结论 如 
下 : 


湿度 条 件 下 ， 随 着 温度 的 升 高 ， 水 分 子 逐 渐 溶 
解 初始 CHO AE, 二 元 C340s+H20 团 簇 相继 出 现 分 
层 结构 (水 分 子 以 团 秘 形式 包 里 初始 团 簇 ) 和 混合 结 
构 (水 分 子 渗入 C0; 团 艇 内 部 ) 。 

不 同 种 类 上 讽 化 钠 盐 的 加 入 对 初始 团 复 的 稳定 性 
(被 水 分 子 溶解 的 难 易 程度 ) 会 产生 不 同 的 影响 ， 其 
稳定 性 排序 为 40 NaF > 40 NaCl > 0 NaX > 40 Nal, 
致使 随 着 温度 的 升 高 包含 Nal 和 NaCl 的 团 艇 出 现 分 
层 结构 和 混合 结构 ， 而 包含 NaF 的 团 艇 只 出 现 分 层 
结构 。 且 这 种 溶解 程度 的 差异 随 着 温度 的 升 高 而 增 
Ke 


高 温 时 ( 产 300 K) ,本 文中 不 同 团 簇 的 表面 均 由 
两 二 酸 分 子 和 水 分 子 构成 ， 且 丙 二 酸 分 子 的 朴 水 基 
朝向 气相 一 侧 。 
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